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Abstract 



^ In this work we consider a simple, approximate, tending toward exact, solution of the 

■ system of two usual Lotka-Volterra differential equations. Given solution is obtained by an 
Tj" , iterative method. In any finite approximation order of this solution, exponents of the corre- 

I sponding Lotka-Volterra variables have simple, time polynomial form. When approximation 

■ order tends to infinity obtained approximate solution converges toward exact solution in some 
, finite time interval. 

o 

"o 

> 

X : 1 Introduction 



As it is well-known [l]-[5] the following system of two nonlinear differential equations of the first 
order ^ 

— = ax — bxy (1) 

^ = -cy + dxy (2) 

represents the remarkable, usual Lotka-Volterra differential equations system with significant ap- 
plications in the ecology, biology, medicine etc.. Here x, y, represent the real, positive variables 
that depends of the time, t, while a, b, c, d, represent real, positive time independent constants, 
i.e. parameters. It is supposed that the initial values 

a;(0) = xo (3) 
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y(0) = ?/o (4) 

of y are known. Existing, exact solution of the system (1), (2) presentable in the simple 
(closed) form (combination of the simple functions of the time, eg. a polynomial of the time etc.) 
is unknown to this day. For this reason different approximate, especially numerical (Runge-Kutta 
etc.), methods for solution of the system (1), (2) are used. But, even in this case approximate 
solution of the system (1), (2) presentable in a simple form is unknown. 

In this work we shall suggest a simple, approximate, tending toward exact, solution of the 
Lotka-Volterra equations system (1), (2). Given solution will be obtained by an iterative method. 
In any finite approximation order of this solution, exponents of the corresponding Lotka-Volterra 
variables have simple, time polynomial form. When approximation order tends to infinity obtained 
approximate solution converges toward exact, simple solution in some finite time interval. All this 
can be very useful for the applications in many domains of the ecology, biology, medicine, etc., 
but explicit consideration of such applications goes over basic intentions of this work. 



2 A simple approximate solution of the usual Lotka-Volterra 
equations system 

As it is well-known [l]-[5] usual Lotka-Volterra equations system (1), (2) can be simply transformed 
in the following system 

^^^^ h fK\ 

a -by (5) 



dt 
d\ny 
dt 

or, after well-known changes of the variables. 



— c -|- dx (6) 



u = lnx X = exp u (7) 

V — Iny y — exp v (8) 
(so that u, V represent the real variables, exponents of x, y), in the system 

^ = a — 6 exp V (9) 



dv 

— — —c + dexpu (10) 
dt 



It is supposed that the initial values 



u{0) = uo (11) 

^;(0) = vo (12) 

of u, V corresponding to initial values xo, yo (3), (4) are known. 

Solution of the system (9), (10) for conditions (11), (12) will be supposed in the following 
approximate form characteristic for n-th approximation step 

u^n) = Ao + Ait + A2t'^ + ... + AJ'' for n = 0,l,2, (13) 

v^n)^ Bo + Bit + B2t^ + ... + 3^^" for n = 0,l,2, (14) 
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where A„ and S„ for n = 0, 1, 2. . . . are unknown real coefficients. Tfiese coefficients will be 
determined by supposition 

Aq = uq (15) 
So = (16) 

and by approximate iterative method 



a — 6exp[t>(„_i)] for n—1,2, (17) 



dt 

— —c + dexp[u(^n-i)] for n—1,2, (18) 



dt 

In the first approximation order, i.e. for n — 1, introduction of (15), (16) in (17), (18) yields 
immediately 

Ai — a — b exp[vo] (19) 

Bi = -c + dexp[uo] (20) 

In this way coefficients Ai and Bi are completely determined. 

In the second approximation order, i.e. for n — 2, introduction of (15), (16), (19), (20) in 
(17), (16) yields 

Ai + 2A2t = a- 6exp[fo + Bit] =a- bexp[vo] exp[Bit] (21) 
Bi + 2B2t = -c + dexp[uo + Ait] = -c + dexp[uo] exp[Ait] (22) 
For additional conditions 

(23) 
(24) 

i.e. for 

(25) 
(26) 

right hands of (21), (22) can be approximated by theirs linear Taylor expansions which yields 

Ai + 2A2t = a-bex-p[vQ]{l + Bit) (27) 

Bi + 2B2t ^-c + d exp[Mo] (1 + Ait) (28) 
It, according to (19), (20), yields 

A2 = -^^I^pN (29) 

B2 = ^^Ii^pN (30) 

In this way coefficients A2 and B2 are completely determined. 

In n + 1-th approximation order introduction of (13), (14) in (17), (18) yields 

Ai + 2A2t + ... + (n + l)A„+ir = a - bexp[vo + Bit + B2t^ + ... + SX] = 

= a-bex-p[vo]ex-p[Bit + B2t'^ + + BJ""] for n = 2,3, .... (31) 



\Bit\ 


< 1 


\Ait\ 


< 1 




1 


i < 








1 


t < 
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El + 2B2t + ... + (n + = -c + (iexp[Mo + Ait + Ast^ + ... + A„r] = 

= -c + (iexp[uo]exp[Ait + ^2^^ + + for n = 2,3,... (32) 

Also, suppose that the following additional approximation conditions are satisfied 

l^\Bit + B2t^ + ... + Bn^l for n = 2,3,... (33) 

1 > > |S2t^| > ... > for n = 2,3,... (34) 

l:^ \Ait + A2t^ + ... + Ant""] for n = 2,3,... (35) 

1 > > l^at^l > ... > for n = 2,3, ... (36) 

Obviously, condition (34) ensures the convergence of (14), while condition (36) ensures the con- 
vergence of the (13). As it is not hard to see, all conditions (33)-(36) can be changed by one 
condition 

t <m^r^(| — 1,1 — — 1,1 — 1,1 — — I) for n-2,3, (37) 

where right hand of the inequality (37) represents the minimum function equivalent to the minimal 
arguments. 

According to (33)- (36), right hands of (31), (32) can be approximated by theirs linear Taylor 
expansions which yields 

Ai + 2A2t+ ...+nAj'' = a-bexp[vo][l + Bit + B2t^ + ... + Bnt""] for n = 2,3, ... (38) 

Bi + 2B2t+ ... + nBnr ^ -c + dex.-p[uQ][l + Ait + A2t'^ + ... + AJ""] for n = 2,3, ... (39) 
Then, it follows 

-bBn-i cxp[vo] r r, o rAn\ 

An — — - for n = 2,3, ... (40) 

n 

dAn-iexp[uo] 

Bn — — - for n = 2,3, ... (41) 

n 



or 



A„= -""'''"°.t"°' -4-2 for n = 2,3.... (42) 



n{n — 1) 

i exp [uq + 
n{n — 1) 



-6rfexp[uo + Wo] , . 

^ri = 7 T\ Bn-2 for n = 2, 3, ... (43) 



Or, 



A2k^ TTzpr, Aq for A; =1,2, 3,... (44) 



m 

= t^^^I^^A, for ^=1,2,3,... (45) 
B,,^tM^^^^±J^B, for ^ = 1,2,3,... (46) 
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(-bdex.p[uo + vo])'' , . 

B2k+i = (2k + iy. k = l,2, 3, ... (47) 

In this way coefficients An and Bn are completely determined by (44)- (47) for n = 2, 3,. . . . 

So, according to previous analysis, coefficients An and Bn are completely determined for any 
n, i.e. for n = 0, 1, 2,.... Obviously, coefficients An and Bn appear explicitly in n-th approximation 
order but their form stand conserved in any next approximation order, for n = 0,1,2,.... 

Further, from (40)- (47) it follows 



Ar. 



A 



n+1 



Bn 



(n + 1) for n = 2,3,... (48) 



(n + 1) for n = 2,3,... (49) 



Bn+l 

It admits that, in a rough approximation, (37) be reduced in 

t^rmn{\^l\^J) (50) 

It demonstrates roughly that suggested approximate solution (13), (14) of the system (9), (10) for 
conditions (11), (12) is consistent and convergent in a time interval [0, r] (whose explicit form will 

not be considered here) that is not infinite small, i.e. infinitesimal. 

Obviously, when n tends toward infinity consistent and convergent approximate solution (13), 
(14) of the system (9), (10) for conditions (11), (12) in a finite time interval [0, r] tends toward 
exact solution. 

Finally, without detailed analysis, we can observe the following. We do not consider here the 
explicit form of [0,r], i.e. r in general case. Also, explicit form of the period, T, of the Lotka- 
Volterra periodical variables x, y is unknown in general case. For this reason we cannot compare 
directly r and T in general case. Nevertheless, we can generally suppose either t > T or t < T. 
In the first case suggested solution (7), (8), (13)-(16), (19), (20), (44)-(47) represents a complete 
solution of the Lotka-Volterra equations system (1), (2) with initial conditions (3), (4). In the 
second case, according to (7), (8), (13)-(16), (19), (20), (44)-(47), we can determine exactly x{t) 
and ?/(t). Then, x{t), y{t) solutions of (1), (2) in a time moment t = r + f for initial conditions 
(3), (4), can be presented by x{t'), y{t') solution of (1), (2) where f' represents new (translated) 
time variable for new initial conditions x{t' = 0) = x{t) and y{t' = 0) = yir). Given solution 
can be obtained by expressions analogous to (7), (8), (13)-(16), (19), ( 20), (44)-(47) and this 
solution is consistent and convergent in some new, finite time interval [0,r']. If r + r' is larger 
than T or equivalent to T obtained solution represents final solution. In opposite case described 
procedure of the solution by "translation" of the time variable must be further repeated. It can 
be supposed that final solution would be obtained after finite number of the repetitions of the 
suggested solution procedure. Simply speaking, exact, simple solution of the usual Lotka-Volterra 
equations would be obtained by a finite number of the suggested analytical continuations. 

3 Conclusion 

In conclusion we can shortly repeat and point out the following. In this work we suggest a sim- 
ple, approximate, tending toward exact, solution of the usual Lotka-Volterra differential equations 
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system. Given solution is obtained by an iterative method. In any finite approximation order of 
this solution, exponents of the corresponding Lotka-Voltcrra variables have simple, time polyno- 
mial form. When approximation order tends to infinity obtained approximate solution converges 
toward exact solution in some finite time interval. All this can be very useful for the applications 
in many domains of the ecology, biology, medicine, etc., but such applications are not considered 
explicitly in this work. 
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